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Abstract 

We construct simple standing wave solutions in a 5D space-time with a ghost-like 
scalar field. The nodes of these standing waves are 'islands' of 4D anti-de-Sitter space- 
time. In the case of increasing (decreasing) warp factor there are a finite (infinite) 
number of nodes and thus a finite (infinite) number of anti-de-Sitter island-universes 
having different gravitational and cosmological constants. This is similar to the land- 
scape models, which postulate a large number of universes with different parameters. 

PACS numbers: ll.25.-w, 11.25.Wx, 11.27.+d 



1 Introduction 

The infinite extra dimension brane theories [H [21 [Sj H] have been one of the most active areas 
of study in the past decade. A concise overview of these models can be found in [5J. The 
original brane models had 5D space-times with 4D (5-function sources which were positive 
or negative tension thin branes. Later works replaced the 5-function sources with localized 
energy momentum tensors of finite thickness which were called thick branes. Examples of 
thick brane models included fluid-like matter sources [6], [7], scalar field matter sources [8],|9], 
and phantom scalar field sources llQ\ [11], [12], |13l HH [15] . Most brane world models focused 
on static solutions, i.e. both metric and any fields being time independent. An exception to 
this are S-brane solutions [lEl HT] HH] [191 [20] which are used for cosmological studies. 

In this paper we present a simple standing gravitational wave solution as a possible 
5D brane model. The existence of this solution requires a negative cosmological constant 
and a phantom/ghost-like scalar field (i.e. a field with a negative sign in front of the field 
kinetic energy term of the Lagrangian) in the bulk. Such ghost-like scalar fields are generally 
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problematic since they tend to make the system unstable. To counter this criticism we will 
show that our model can be embedded in a 5-dimensional Weyl model of gravity [2T1 ESI 
[21]. Because of the geometrical origin of the coupling of the scalar to the space-time and 
since the Weyl scalar does not couple with other matter fields we are able to avoid the usual 
problems of instability due to ghost fields. 

We assume that the ghost-like field vanishes on the brane at the origin. This condition 
leads to a quantization of the oscillation frequency of the standing wave which in turn makes 
the scalar field vanish at a finite or infinite number of nodes of the standing wave depending 
on whether the warp factor of the metric is increasing or decreasing. By taking the matter 
fields to bind to the standing waves nodes one has a finite or infinite number of anti-de- 
Sitter island-universes, each with different parameters. In the simple model presented here 
the different parameters are the effective 4D gravitational constant and the 4D cosmological 
constant. These two 4D constants are fixed in terms of the 5D constants and the value of 
the warp factor at the nodes. 

For the case of an increasing warp factor, and a finite number of anti-de-Sitter 'islands', 
one could address the family puzzle by taking different fermion generations as bound to dif- 
ferent nodes and fixing the number of nodes in the bulk at three. For the case of a decreasing 
warp factor, and infinite number of anti-de-Sitter island-universes, one has something like the 
landscape picture in string theory/M theory [251 [2Sl 123) with its large number of universes 
with different parameters. 

In the model presented in this paper ordinary matter fields are assumed to be bound 
to one of the anti-de-Sitter 'islands'. This is an advantage over the original one-brane or 
two-brane models of [H |21 El S] , where the ordinary matter fields are assumed to be localized 
to a brane with positive or negative tension. In these models one needs to explain why these 
brane tensions are not observed. 

There are several open questions or problems with the present 5D standing wave solution 
as brane world model: The first question concerns the stability of the whole 5D solution 
given the presence of the ghost-like scalar field. This issues of the stability of brane world 
models in the presence of sources which violate some or all of the usual energy conditions 
(e.g. negative tension branes [U 121 [3l H] , or ghost fields [lOl [TTl [T21 [13] ) has been known since 
the beginning of the study of brane world models [28]. In section 3 we present a possible 
resolution to this issue by considering Weyl's generalization of Riemannian geometry. We 
also offer some discussion in the conclusion section of cases where a ghost field, at least at 
the classical level, leads to greater stability of brane world models as compared to normal 
scalar fields. The second question is: "On or near the nodes are the gravitational fluctuations 
confined so as to give effective 4D gravity near the nodes?" . In section 6 will show that near 
the nodes one does get effective 4D gravity in much the same way as in [H [21 [1] . The last 
question/problem is give a localization mechanism for the matter fields (i.e. fields with spin 
0, spin |, spin 1) on the brane. This localization problem has been an unresolved question 
for brane world models in general [23 [301 Ell E21 E3l IMl |35]. Here we simply postulate that 
matter fields are bound to the nodes of the standing wave solutions. In section 7 we give two 
possible mechanisms for accomplishing this binding of matter fields to the nodes. The first 
mechanism involves the coupling of the quadrupole moment of the stress-energy tensor of 
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the matter fields to the Riemann tensor. The second mechanism involves coupling of matter 
fields to the ghost-like field. 



2 Metric and Einstein's equations 

We study solutions of the 5-dimensional Einstein equations, 
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Rab — T^QabR — — ^ — Tab — J^bQab , (1) 
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where G5 and A5 are 5D Newton and cosmological constants respectively [11 |2l |3l H] . Lower 
case Latin indices a, h run over the full 5D space-time 0, 1, 2, 3, 4. The metric ansatz we use 
is 

ds^ = e'"!''! {dt^ - e"dx2 - e'^dy^ - e-^'^dz^) - dr^ , (2) 

where a is a constant and the ansatz function u = u{t,r) only depends on the time t and 
the extra coordinate r. This ansatz is a combination of the 5D warped brane world model 
through the e^"''"' term [H |2l |3l H] plus an anisotropic {t, r)-dependent warping of the brane 
coordinates, x, y, z, through the terms e"^*'''-', e"^"^*'''-'. The absolute value sign around r gives 
a 5-function source/brane at r = exactly like the one brane models of [H 121 El S] • Since we 
will focus on standing wave solutions 'to the right', i.e. r > 0, we drop the absolute value 
sign in ([2|). We study both decreasing (a < 0) and increasing (a > 0) warp factors. We 
also find that for both increasing and decreasing warp factors the ansatz function u{t, r) has 
nodes - u{t,r) = - at specific values of r. At these nodes the space-time in ([2|) reduces 
to 4D Minkowski space-time plus a scaled, negative cosmological constant - at the nodes 
the space-time becomes effectively 4D anti-de-Sitter. As one moves away from these nodes 
there is an anisotropic stretching/shrinking along the x,y,z directions due to the metric 
components e"^*''"-', e~'^'^^^'^\ If u{t, r) is positive (negative) as one moves away from the node 
in r the x, y directions will expand (shrink) as while the z direction will shrink (expand) 
as e~^". In section 7 we will use this feature of the metric to propose a possible localization 
mechanism for matter fields. 

For the ansatz (|2|) we find the non-zero components of the Ricci tensor: 

2ar+u f ig-2ar^ _ 2a^i' - -u" - Aa^ 
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2 2 
{-e-^^'u + Aau' + u" -Aa^) , (3) 

-4a2 , 



where overdots mean time derivatives and primes stand for derivatives with respect to r. In 
([3|) there should be terms proportional to 6{r) coming from the tension of the thin brane at 
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r = 0. We do not write this explicitly since we focus on r > and have thus dropped the 
absolute value in the warp factor of the metric ([2]). 

3 Source scalar field 

To the metric above we add a massless, non-self interacting scalar phantom/ghost-like field, 
(pit,z) which obeys the Klein-Gordon equation, 

^ daiVag^'M) = e-"''4> - <\>" - 4a0' = . (4) 
V 9 

Here g is the determinant of the 5-dimensional background space-time given by ([2]). The 
energy-momentum tensor of the phantom-like field is taken in the form: 

Tab = -da<i)db(t) + ]^gabd''(i) 5c0 • (5) 

Strictly speaking cj) is not a phantom field as defined in [38] , where the criterion for a phantom 
field was p/p < —1 {p and p are the pressure and energy density of the field respectively). 
From (IS]) one can obtain p and p for the field and since the field is non-self interacting one 
does not have pj p < —1. 

To avoid the well-known problems with stability that occur with ghost fields we can 
associate the ghost-like field with the geometrical scalar field in a five- dimensional, inte- 
grable Weyl model. In Weyl's model a massless scalar appears through the definition of the 
covariant derivative of the metric tensor, 

Dcgab = gabdc(t) ■ (6) 

This is a generalization of the Riemannian case where the covariant derivative of the metric 
is zero. The result in ([6]) implies that the length of a vector is altered by parallel transport. 
Weyl's scalar field in (|6]) may imitate a massless scalar field - either an ordinary scalar or 
ghost-like scalar [211 [22l [23|, [21] • The gravitational action for Weyl's 5D integrable model 
can be written as 

= TTtV / ^'^v^ - (6 - 50 ^'^^ 9,0] , (7) 
ib7rG5 J 

where ^ is an arbitrary constant which can take any value. For example, if one takes 
^ = 13/10 in ([7]) then the coefficient in front of (9"0 9,0 becomes —1/2. This would give a 
ghost field which would lead exactly to the equation of motion (jlj) and the energy-momentum 
tensor ([5]) for our phantom-like scalar field. Thus we can start with a 5D Weyl model and 
require that we have Riemann geometry on the brane by assuming that the geometrical 
scalar is independent of brane coordinates and vanishes on the brane. We will show 
later that our solution has exactly this character - the scalar field vanishes at the nodes of 
u(t,r) and is independent of x,y,z. Associating our scalar field with the geometrical Weyl 
scalar defined via (E]) avoids the usual instability problems of ghost fields since the Weyl 
model is known to be stable for any value of ^. 
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4 The solution 



For ([5]) one can rewrite the Einstein equations ([T]) in the form: 

3 

Rab = -datrdbO- + -QabK , (8) 

where the gravitational constant has been absorbed via the redefinition of the scalar field, 

a = ■ (9) 

By combining ([8]) and ([3]) one can see that the "constant" terms from i?^,^ (i.e. those terms 
that up to some metric factor are ±4a^) can be canceled if the 5D cosmological constant is 
chosen as 

A5 = . (10) 

This is the same as the fine tuning used in standard 5D brane models [H [21 [3l H] . 

The terms from Rf^,y which are quadratic in u can be accounted for if one assumes that 

^it,r) = ^uit,r) , (11) 

so that u{t,r) satisfies (jl]). A similar equality between the scalar field and metric ansatz 
function was found for the domain wall plus standing wave solutions of [SH] • Since our scalar 
field is proportional to the metric ansatz function u{t, r) the scalar field will vanish wherever 
u{t,r) has nodes. This requirement, that the scalar field vanish at the zeros of u{t,r), was 
one of the necessary conditions pointed out in the last section in order to be able to associate 
our scalar field with a Weyl scalar field. 

Because of ( ITOl) and ( |TT|) . solving Einstein equations has been reduced to finding solutions 
to the ordinary differential equation, 

u-u"- Aau' = . (12) 

We want to have a standing wave solution to (IT^ . so we use separation of variables 
writing 

u{t,r) = Csm{ujt)f{r) , (13) 
where C and u are constants. Because of flTT]) the same separation applies to cr, but with a 
different constant, C — t- \/^C. Equation ( !T2|) now becomes ( !T3|) is: 



/" + 4a/' + w2e-2«7 = . (14) 
The general solution to this equation is: 

/(r) = J2 (^e-*^") + Be-'^'N^ (^e""") , (15) 
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where A, B are constants and J2 and N2 are 2 order Bessel functions of the first and second 
kind respectively. Normally the N2 Bessel functions are discarded since they blow up at the 
origin. But here the functional dependence is e~"^, rather than r, and neither J2 nor A^2 
diverges at r = 0. 

Before moving on to the discussion of the physical meaning of the solutions in ( !T5|) we 
add the boundary condition that the ghost-like field should vanish at the brane, r = 0. To 
accomplish this we should take either A = or B = 0, since the zeros of J2 and N2 do not 
coincide. Then we set 

- = X2,n , (16) 

a 

where X2,n is the n'^^ zero of the 2"°' order Bessel function J2 or N2, depending on whether 
one takes A = or 5 = in (JT5i) . The condition ( |T6i) quantizes the oscillation frequency, uj. 



5 Standing waves as Anti-de-Sitter 'islands' 

We now analyze the physical consequences for metric given the solutions given by (IT^ . 
and ffT^ . 

First we note that ([2]) in terms of the behavior of the x,y,z coordinates gives different 
asymptotic properties for the increasing (a > 0) and decreasing (a < 0) warp factor cases. 
For the increasing warp factor (a > 0) the metric function, u{t,r) oc f{r), decreases like 
e"^'"' as one moves into the bulk. As one goes to large r — > 00 this factor of e"^"'' drives 
u{t, r) — )■ and one has an anti-de-Sitter space-time but scaled up by an overall factor e^"^. 
For the decreasing factor (a < 0) the function u{t, r) grows like e^""^ into the bulk therefore 
distances in the x and y directions expand like Exp(e^"''), while in the z direction they shrink 
like Exp(-2e2"''). 

The second observation is that u{t, r) is oscillatory and has zeros at various points along 
the r-direction - whenever /(r) = 0. This happens when J2{f') = or N2{r) = 0, depending 
on if one is considering B = or A = 0. At these zero points, r^, one has u{t, r^) = and 
the metric in ([2]) reduces to the standard 5D brane metric with an exponential warp factor 
[H 121 m H]. This spatial oscillatory behavior of u{t,r) leads to 'islands' of anti-de-Sitter 
space-times in the bulk. 

The third observation is that the physical parameters on each node are scaled by the 
warp factor e"'' or some power thereof. For example, taking the 4D reduction of 5D Einstein 
equations ([1]) one finds that on the nodes of the anti-de-Sitter 'islands' there are effective 
4D, negative cosmological constants given by 

A4 = e^'^'^'-As . (17) 

Thus the values of the effective 4D cosmological constant on the nodes is set by the scaled 
5D cosmological constant. In the usual thin brane models [H [21 |3l H] the 5D cosmological 
constant is fine tuned to exactly cancel the 4D brane tension. Thus on the brane the effective 
4D cosmological constant is zero. In our case at the nodes, r^, there is no brane tension and 
so our effective 4D cosmological constant is non-zero and given by (fTTj) . For a decreasing 
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wrap factor (i.e. a < 0) one has a cosmological constant which is exponentially suppressed 
relative to the true 5D value, A5. Below we argue that these nodes are 4D anti-de-Sitter 
island universes on which matter fields are bound and on which gravity is effectively 4D. 

For the case a > and B = 0, there will be (n + 1) such 'islands' (n is the number of 
the zero from (fT6|) ). As r runs from to 00 the argument of J2 in (fT5|) runs from X2,n to 
giving n zeros. One additional zero comes from /(r — )■ 00) = 0. The zeros will occur at 
Tm where u:e~°'^'^ /a = X2^m with X2^rn being a zero of J2 with m < n. The case a > and 
A = works out the same with (n + 1) 'islands' since one also has an additional zero coming 
from /(r — 7- 00) — ?■ - in this case the divergence coming from A'"2(0) is dominated by the 
g-2ar pre-factor in f|T5l) . The ue~"'^ /a dependence of /(r) has the effect of stretching out the 
zeros of the Bessel functions, i.e the spatial oscillation frequency decreases with r. 

In the case a < (with either A = or B = 0) there are an infinite number of zeros 
for u{t,r) due to the ue°''^/a dependence of f{r). In addition the zeros are compressed as r 
increases, i.e. the spatial oscillation frequency increases with r. This compression (a > 0) 
and stretching (a < 0) of the location of the nodes has a connection with the effective 
thickness of the 4D anti-de-Sitter islands - the compression of the nodes tends to lead to 
'islands' with a smaller effective thickness while stretching of the nodes leads to 'islands' with 
a large effective thickness. We discuss this further in section 7 where we examine localization 
mechanisms for matter fields to the nodes. 

6 Localization of gravity 

Near r = in metric ([2]) one has a node+brane as in the 5D single brane models of [H |2l |3l H] . 
The gravitational perturbations around r = will have one delta-function bound state and 
continuum states which start from zero mass i.e. which do not have a mass gap. Thus at 
the r = node one has effective 4D gravity as in the original models. We now look at the 
nodes of u{r,t) for r > to see if near these nodes gravity also becomes approximately 4D. 
This is done by studying the 4D perturbations of the metric near the nodes. If there is a 
zero mass graviton mode near the node then its exchange between particles localized on the 
node will lead to a Newtonian potential with corrections coming from the massive modes 
(m > 0). Close to the nodes one has u{r,t) ^ and the metric ([2]) takes on the usual 5D 
warped geometry. Near the nodes small fiuctuations around this background can be written 
as 



where = (t,x,y,z) and Greek indices run over 4D space-time - /i, 1/ = 0,1,2,3. The 
tensor, /i^,^, gives 4D perturbations around the usual 5D brane world background. Further 
we fix the gauge so that h^^, is transverse and traceless. 




(18) 



(19) 



Next we separate the tensor perturbation into a ID and 4D part as 



(20) 
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with /i|tt^(x^) satisfying 

=p' = m2 . (21) 
Putting all this together yields the following equation for the ID part of the perturbations 

^" - Aa^tP - m^e-^^^ij = . (22) 

In order for gravity to be effectively 4D near the nodes fl22l) should have a zero mode solution 
- m = 0. It is easy to see that it does have a zero mode given by, 

Mr) = e"'^ , (23) 

up to a normalization constant. Equation also has m ^ solutions which are similar to 
( |T5|) i.e. combinations of second order Bessel functions of the first and second kind, J2,N2, 
[5]. The exchange of the zero mode, m = 0, between massive particles fixed on the node, 
will lead to a 1/r Newtonian potential while the exchange of the m ^ modes will lead to 
corrections which go as l/(a^r^) [H |2l [Sj HJ |5] . 

The effective 4D Newton's constant on a particular node can be obtained using the above 
results. Let us focus on one particular node, r^^. The massless gravitational perturbation is 

h,. = e'-^-h';f}{xn , (24) 

where h^^}{x^) is the 4D part. The effective 4D Newton's constant can be obtained by 
inserting fl2^ into the 5D gravitational action and looking at the quadratic part 

Sg = -i— r^'±_ [ d^^^dhf 



1 



16nG. 



Tm+d 



~d 



dr e^'^' / <fx{dh^'^^f (25) 



where 2d is the width around the node, r^, within which the matter particles are assumed 
to be bound. The index m satisfies m < n where for m = n one is located at the r = 
brane. Using the last expression in (125|) the effective 4D Newton's constant can be written 
as 

G — '^^5 ^ G5 

The last expression assumes a small thickness i.e. c? ^ 1, or ad <^ 1. In the case when a > 
one can exponentially suppress the 4D Newton's constant, G4, relative to the 5D Newton's 
constant, G5. 
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7 Localization of matter 



So far in this paper we have assumed some mechanism for binding matter fields to the nodes 
of the standing waves. Now we want to give two possible localization mechanisms. 

The first mechanism is borrowed from condensed matter physics. It is known that stand- 
ing electromagnetic waves, so called optical lattices, can provide trapping of various particles 
by scattering and dipole forces [101 SI] . In 113] localization was demonstrated through 
quadrupole forces as well. It is also known that the motion of test particles in the field of a 
gravitational wave is similar to the motion of charged particles in the field of an electromag- 
netic wave IH] . Thus standing gravitational waves could also provide confinement of matter 
via quadrupole forces. As an example let us consider the equations of motion of the system 
of spinless particles in the quadrupole approximation HHl HZ] , 

^ = ^^w = —^r'''D^R^,,s , (27) 

where is the total momentum of the matter field/particle and J°^'>'^ is the quadrupole 
moment of the stress-energy tensor for the matter field/particle. The oscillating metric due to 
gravitational waves should induce a quadrupole moment in the matter fields. If the induced 
quadrupole moment is out of phase with the gravitational wave the system energy increases 
in comparison with the resonant case and the fields/particles will feel a quadrupole force, 
^quadi which cjects them out of the high curvature region i.e. it would localize them at the 
nodes. This gravitational binding mechanism would be effective for all types of matter fields 
since gravity couples to all forms of energy-momentum. Since the matter fields are ejected 
from the high curvature region and toward the nodes the width of the anti-de-Sitter 'islands' 
(i.e. the thickness of the branes) depends on how rapidly the ansatz function u{t,r) oc /(r) 
from (fT3l) ( |T5l) changes from zero. The width of some anti-de-Sitter 'island' will depend in 
the derivative of /(r) at the node - /'(r^). The larger /'(r^) is the smaller the width of the 
anti-de-Sitter 'island'. Now /'(r^) depends on the size of either A or B and as well as a, u 
and Tm- For some choice of these parameters some nodes might not be phenomenologically 
acceptable since the brane thickness might be too large. Also in general for the case when 
a > there will be fewer phenomenologically acceptable branes since in this situation the 
nodes are stretched out, as discussed in section 5, and this tends to decrease fir^) as r 
increases. On the other hand when a < there is an increase of f'{rm) as r increases which 
leads to more nodes which have a phenomenologically acceptable thickness. 

The second mechanism would be to propose some coupling between the ghost-like field, 
(T, and matter fields of the form aAA, where A is a scalar, spinor, vector or tensor field. For 
a normal scalar field such a coupling leads to an attractive force, while for a ghost-like field 
it leads to repulsion of the matter fields from regions with non-zero a. This would force the 
fields A to regions were a vanishes, i.e. the nodes of the standing waves. Note 

that coupling ordinary matter fields to ghost fields is problematic since in general it leads 
to instabilities. However, as mentioned in the Sec. 3, our massless ghost-like field can be 
associated with the geometrical scalar of integrable Weyl models. 

Both binding mechanisms have open questions (e.g. "Does the gravitational standing wave 
induce a quadrupole moment of the matter field stress-energy tensor and if so what is its 
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form?"; "Can one safely couple ordinary matter fields to ghost-like fields?"). A detail an- 
alyzes of localization of specific matter fields via these two possible mechanisms will be 
considered in a forthcoming paper. 

8 Conclusions and discussion 

A simple standing wave solution for 5D spacetime with a ghost-like scalar field plus 5D 
negative cosmological constant was presented. The requirement that the ghost-like field 
vanish on the brane quantized the standing wave oscillation frequency in terms of the zeros 
of the 2"''^ order Bessel functions. Thus the ghost-like field is not observable on the brane, at 
r = 0, nor at any of the nodes, r^, of the standing wave. This is similar to the static brane 
model [ini HI 112], where the phantom/ghost field only existed in the bulk. 

For the case of increasing (decreasing) warp factors the ghost-like field, a(t,r), and the 
metric function, u{t,r), vanish at a finite (infinite) number of places in the bulk forming 
'islands' of 4D anti-de-Sitter space-time. By assuming that the matter fields are bound to 
these nodes we arrive at a model where each node is an anti-de-Sitter island-universe with 
different parameters, i.e. an effective 4D Newton's constant from fl2B]) and an effective 4D 
cosmological constant from f[T7l) scaled by the value e^"^'™ at the particular node. One could 
address the hierarchy problem and the small size of the cosmological constant by taking 
our universe as one of these Bessel nodes, where the 5D Newton and cosmological constants 
are appropriately scaled to their effective 4D values. A problem is that the scalings go in 
oppose directions. For example, with a > (1171) gives a 4D cosmological constant which 
is exponentially larger than the 5D cosmological constant, while (126 p gives a 4D Newton's 
constant which is exponentially smaller than the 5D Newton's constant. For a < the 
effective 4D Newton's constant and the effective 4D cosmological constant have opposite 
scaling to the a > case. 

There are two distinctly different cases for the standing wave solutions - increasing warp 
factor (a > 0) and decreasing warp factor (a < 0) - corresponding to a finite, or infinite num- 
ber of anti-de-Sitter 'islands' respectively. The finite case could be applied to the generation 
problem by fixing the model to have three 'islands' with different fermion generations bound 
to different nodes. This is different from brane world models of the generation puzzle like 
[IHlllHlEn], where the generations were associated with different zero modes bound to a single 
brane. For the infinite branes case one has a simple version of the landscape picture coming 
from string/M theory, where one has a large number of anti-de-Sitter island-universes each 
having different parameters e.g. effective 4D gravitational and 4D cosmological constants. 

The present models has the same kind of warped geometry as the usual brane models 
[H [21 [3111]. In distinction from these models the "branes" (i.e. the nodes of M(r, t), excluding 
the one at r = 0) of the present standing wave background solution do not have a brane 
tension but are simply 4D anti-de-Sitter space-times. This is an advantage since for the 
models of [H [21 [3l [1] one should explain why the (5-source brane tension is not observed. 
The disadvantage of the present model is the need for a ghost-like field. The need for the 
ghost-like field is ameliorated by the fact that it vanishes on the r = brane and on all the 
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Bessel node anti-de-Sitter 'islands'. 

The issue of instability due to an unusual matter source is also a problem for the original 
two-brane model which has a negative tension brane. If the negative tension brane is free to 
oscillate this results in arbitrarily large negative energy modes making the system unstable 
[28] . A similar problem occurs in the present model if we take our scalar field as an ordinary 
ghost field. However, if we associate our ghost field with the geometrical scalar of a 5D 
integrable Weyl model [221 1231 121] this alleviates some of the instability problems since 
the scalar field comes from the metric via In addition for Weyl models the ghost-like 
field does not have the standard couplings with the brane energy-momentum and thus does 
not cause instabilities. 

Additionally previous work with standard ghost fields shows - that at least at the classical 
level - bulk ghost fields may in fact be better at stabilizing brane world models as compared 
to bulk regular scalar fields. In [^I] bulk ghost fields in 5D are shown to lead to radion 
stabilization for models with positive tension branes. Essentially the bulk ghost field replaces 
the negative tension brane. Further it was shown that a bulk scalar field whose sign can 
vary between regular (positive) and ghost (negative) can be used to model both positive and 
negative tension branes. Also, there are certain configurations of the bulk ghost fields of 
reference [51] which lead to stronger localization of gravity than the original brane models. 
In references [521 IS3] it is shown that the two brane model of [U [21 El II] is stabilized by a 
ghost /tachyon bulk field, while for a regular bulk scalar field the two brane model is unstable. 
While the stability issue is crucial for all brane world models it has not yet found a complete 
solution even for the original models of [H [21 [3l [1]. Here, as in [TOl [HI [T21 [51j we simply 
use the ghost scalar to build a brane world model from a 5D standing wave, leaving the 
generally unresolved question of stability of the solution for a future work. As a final note 
in [5^ [551 [56] it is shown how to construct a consistent field theory with tachyons in the 
context of D-branes. 

Although throughout this paper we have assumed some mechanism for binding mat- 
ter fields/particles to the anti-de-Sitter 'islands' in section 7 we put forward two possible 
mechanisms for accomplishing this binding. The first localization mechanism was based on 
the quadrupole force in ([27]) which would eject particles from the high curvature regions 
and thus drive them to the nodes. One might give a heuristic description of this binding 
mechanism by noting that as one moves away from the nodes where u{t, r) = that the 
brane coordinates x, y, z are distorted in an anisotropic way - the x, y coordinates will be 
stretched/compressed while the z coordinate will be compressed/stretched. This anisotropy 
causes a force, via ( 127]) . which tends to drive the matter fields toward the nodes. The sec- 
ond mechanism also involved ejecting the matter fields from the anti-nodes and towards the 
nodes, but in this case the ejection mechanism was accomplished by coupling the matter 
fields to the ghost-like scalar. This would give rise to a repulsive force which would drive the 
matter fields toward regions where the ghost-like scalar field was small i.e. toward the nodes. 
For both mechanisms the thickness of the brane was related to how rapidly the metric and 
scalar field ansatz functions changed near the nodes i.e. to have a thinner brane one should 
make firm) larger. In general the case with a < would have thinner branes. 
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